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Abstract 

We study the collapse of two-dimensional polymers, via an 0(n) model on the square lattice that 
allows for dilution, bending rigidity and short-range monomer attractions. This model contains two 
candidates for the theta point, 0bn and 0ds, both exactly solvable. The relative stability of these 
points, and the question of which one describes the ‘generic’ theta point, have been the source of a long¬ 
standing debate. Moreover, the analytically predicted exponents of 0bn have never been convincingly 
observed in numerical simulations. 

In the present paper, we shed a new light on this confusing situation. We show in particular that 
the continuum limit of 0bn is an unusual conformal held theory, made in fact of a simple dense polymer 
decorated with non-compact degrees of freedom. This implies in particular that the critical exponents take 
continuous rather than discrete values, and that corrections to scaling lead to an unusual integral form. 
Furthermore, discrete states may emerge from the continuum, but the latter are only normalizable—and 
hence observable—for appropriate values of the model’s parameters. We check these hndings numerically. 
We also probe the non-compact degrees of freedom in various ways, and establish that they are related 
to huctuations of the density of monomers. Finally, we construct a held theoretic model of the vicinity 
of 0BN and examine the how along the multicritical line between 0bn and 0ds. 


1 Introduction 


Two-dimensional self-avoiding walks with additional generic short-range attraction are believed to experience 
a collapse transition as the temperature is lowered |^[^. The high-temperature phase is in the universality 
class of ordinary self-avoiding walks (SAWs), also known as dilute polymers^ with well-known exponents. 
The low-temperature phase is in the universality class of so-called dense polymers, also with well-known 
exponents IQ In between sits the so-called theta point, where the polymers are somewhat more compact 
than in the dilute phase, but still have a fractal dimension strictly smaller than two. 

The theta point is well identified in the language of the 0(n) held theory (in the polymer limit n —>■ 0) 
as a tricritical point 01 . while the ordinary SAWs correspond to the usual critical point. Non-generic 
attraction between monomers would lead to higher criticality for the transition point, leading to what is 
sometimes called a theta prime point 


10,11 


It is unfortunately difficult to identify ‘the’ theta point based on its expected tricritical nature. When 
dealing with geometrical problems, the counting of physical observables—many of which are, in a sense, non¬ 
local— is ambiguous. Moreover, in two dimensions, the Landau-Ginzburg picture is hardly able to organise 
the zoo of known universality classes, which is rendered even more complicated by the lack of unitarity 


^The zero-temperature case, where the walks are fully packed, is known as Hamiltonian walks. It presents more subtle 
features and strong non-universality, with critical exponents that are lattice-dependent and depend continuously on the 

stiffness for models of semi-flexible walks [^|^. We will not discuss this further here. 
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inherent to geometrical problems. Finally, producing reliable numerical results turns out unexpectedly hard. 
This means the identification of the theta point critical exponents has led to a long-standing controversy. 

were able to solve exactly a model (first proposed by Coniglio et 


In 1987, Duplantier and Saleur 12 


al. [^) of SAWs on the honeycomb lattice with a particular type of attractive interaction. They conjectured 
that this interaction was generic enough to put the model in the theta universality class, and obtained the 
exponents 


ve 

7e 


0.571, 

1.143, 


7 


0.429. 


( 1 ) 


Later, the same authors [14| carried out extensive numerical simulations to explore the stability of their 
conjectured theta point against additional attractions, and concluded that indeed, the exponents they had 
obtained described the generic (tricritical) theta point. These exponents were also found in agreement with 
numerical simulations of interacting self-avoiding walks on other lattices 15 , but it is fair to say that the 


proximity of the dense phase with very strong corrections to scaling rendered these results a little less definite 
than one would have liked. For the square lattice and presumably generic nearest-neighbour attractions, the 
best recent numerical results are 


16 


num 


= 0.570(2), 
= 0.46(3). 


( 2 ) 


Meanwhile, in a beautiful series of papers 17 -20 , Nienhuis and collaborators managed to find an exactly 
solvable 0(n) model on the square lattice with a rich set of phase transitions. One of these has, for u —>■ 0, 
all the required characteristics of the theta point, but leads, disturbingly, to critical exponent^ which are 
different from those of 12 


12 

t'BN = ^ ~ 0.522 , 

53 

7BN = 1-152. 


(3) 


It was suggested at the time in 20 that this solvable point might be in fact ‘the’ generic theta point, implying 
that the point obtained by Duplantier and Saleur 12 was, in fact, a higher multicritical point - although it 
is hard to see how the exponent z/ in Q for a higher multicritical point could be larger than the one ([^ for 
the generic theta point. 

To complicate the matters further, the exponents obtained in 


20 , despite the impeccable derivation 


obtained via the Bethe Ansatz, have never quite been observed in numerical simulations. Frustratingly, for 
the very model discussed in 20 , and using the most state-of-the-art numerical techniques available today, 
the exponents obtained in direct simulations 22 read 


■gr = 0.576(6), 
= 1.045(5). 


num 

7bn 


(4) 


and are, bizarrely, close to the ones Q obtained by Duplantier Saleur! Similar conclusions were drawn 
21 in the context of the measure of surface critical exponents, both from transfer matrix and DMRG 


calculations. 

The discrepancy with respect to the analytical results (§ is essentially unheard of in the field of exactly 
solvable models and conformal field theory in two dimensions, where, usually, theory and numerical exper¬ 
iments match almost perfectly]^ Of course, one important point is that in 20 , the quantities which are 


^The subscript BN refers to Blote and Nienhuis |17| and will be used extensively throughout this work. 

®To give but one example of an (at least seemingly) very comparable situation, the exponent 7 for Hamiltonian walks on the 
square lattice was found numerically to be 7 = 1.0444(1), in perfect agreement with the subsequent analytical solution [q 

7 = Hi ~ 1-0446. 
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Figure 1: Configuration of the semi-flexible vertex-interacting self-avoiding walk (VISAW). 


studied involve a grand canonical ensemble of walks where the length fluctuates, while in the simulations 
of 22 polymers have a fixed length, and the ensemble is canonical. However, this point is usually easily 


taken care of by considerations of Legendre transform, and has never, in all the other cases studied so far, 
led to any particular difficulty. 

While evidence gathered over the last many years (see, e.g., 23p4 for recent contributions) has confirmed 


that the exponents obtained in 
solved model in 


20 


12 


indeed are those of the generic theta point, the meaning of the exactly 
has remained totally unclear: what is its universality class in the 0(n) field theory 
language? Why are the exponents so hard to observe numerically? 

The detailed answer to this question is the object of this paper and turns out to be a rather long story. 
It can however be summarised rather succinctly: The universality class of the model in 20 is profoundly 


different from the one of ordinary SAWs. While for the latter, the spectrum of critical exponents is discrete — 
as is the case for most familiar model^of statistical mechanics—for the former, the critical exponents form 
in fact a continuum, of which the values (|^ obtained in [20| are simply the lower bounds. While, as we shall 
see below, this explains the numerical difficulties encountered so far, it also shows that the continuum field 
theoretic description is profoundly different from the 0{n —)■ 0) model theory expected for the theta point. 
The field theoretic description instead involves a non-compact target, and is related with the Black Hole 


sigma model conformal field theory that has been discussed intensively in the string theory literature 30 -32 . 


2 The polymer model and its phase diagram 


The model considered in is one of a single polymer on the square lattice, where edges can be visited at 
most once. The polymer is not allowed to cross itself, but collisions where two pieces of the polymer barely 
avoid each other are allowed at all sites: an energy —et is associated to this in order to take into account the 
presence of medium-range attractions in the physics of the theta point (see figure . 

Additionally, it is convenient to also allow for some stiffness, and thus associate with two parallel con¬ 
secutive monomers an energy —eg. The resulting object is called a semi-flexible VISAW (vertex-interacting 
self-avoiding walk), in the notations of 


22 that we follow here. The partition function of a polymer made 


of N monomers is thus given by 


Zn{t,p) = 


VISAW 


^number of doubly visited sites ^number of straight segments 


(5) 


with T = and p = . While Monte Carlo simulations focus on polymers of fixed large length N, the 

integrable model deals instead with a grand canonical ensemble where the monomers have a fugacity K, so 
the partition function is 


OO 

N^O 


N 


Zn(t,p) . 


( 6 ) 


The phase diagram in the {p, K, r) parameter space, which has been the object of several numerical 
studies 33 -35 , is depicted schematically in hgure For small values of the interaction parameter t, the 


few noticeable exceptions have nevertheless appeared in the recent literature phteQl. 
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Figure 2: Schematic phase diagram of the two-dimensional VISAW. The nature of the different phases and 
transitions between them is explained in the main text. The two candidates, 0ds and ©bNi for the theta 
point universality class are represented, respectively, by an orange (DS) and a green (BN) dot. 


physics is essentially the same as for usual (non-interacting) SAWs: for small values of the fugacity K the 
polymer is in a massive phase and remains of finite total length as the size of the system is sent to infinity, 
whereas for high values of K the polymer is in a critical, dense phase where it covers a finite fraction of the 
lattice edges in the thermodynamic limit; the two phases are separated at intermediate AT by a different, 
dilute critical phase (represented by a red surface in figure]^. 

In order to investigate the effect of larger values of r, it is instructive to consider Ising degrees of freedom 
living on the faces of the square lattice, such that two Ising spins on adjacent faces share the same orientation 
iff they are separated by a monomer. While this Ising degree of freedom is disordered in the usual SAW 
dense phase, it is clearly ordered in the so-called completely packed (Eulerian walk) limit of large K and r, 
and therefore some critical Ising surface should be expected in the phase diagram (shown in blue in figure 
[^. Above this critical surface, the transition between the massive and dense (ordered) phases is now of first 
order, and is represented by a dashed surface in the figure. 

Whereas all the phases described so far are by now well understoodj^this is not the case of the multicritical 
line joining the dilute, Ising critical, and first-order surfaces, and which contains both candidates under 
consideration for the theta point physics: 


• The integrable point ©bn 
throughout this paper. 


20 , depicted as a green dot in figure which will be our main focus 


The universality class of 12 is known to describe the particular integrable point ©ds situated at 
p = 0, that we depict as an orange dot in figure 


The precise location in the (p, AT, r) phase diagram of the points ©bn and ©dSj as well as that of three 
further integrable points, follows from the integrable solution |17| and various transformations, of which 


we shall review the necessary details in section Suffice it here to say that the integrable model can be 
formulated in terms of a loop model with the following nine different local configurations around each lattice 
site: 


(7) 



®Note however that the p = 0 plane reveals some surprises for generic n [35] that we plan to unravel in a future piece of 
work |36|. 
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The corresponding Boltzmann weights are pi,..., pg as indicated, and there is an additional non-local weight 
n per loop. We set the loop weight n = 0 to obtain polymers. The integrable model gives the same weights to 
local configurations that are related by a horizontal and/or vertical reflection, so that pg = Pa and p 4 = ps. It 
contains a further (so-called spectral) parameter and various gauge factors ±1 that can be adjusted to obtain 
isotropic solutions for which we have also pg = P 4 and pg = pg. The correspondence between (pi, pg, pe, ps) 
with the parameters used in hgurej^is then easily seen to be 


Pg 

P = — , 
P2 


K=^ 

Pi 


T = 


PlPS 

{P2? 


( 8 ) 


We also note that 22 prefers to trade r for another parameter defined by w = K^r = ps/pi- 


The integrable model studied in 17-^ is defined in terms of trigonometric weights (pi,... ,pg;n) whose 
periodicities are such that, after constraining to isotropic n = 0 cases, there are in fact four distinct solutions 
for the parameters (p,K,t). These are referred to as branches 1, 2, 3, 4 in [^, and as regimes I, II, III 
in (regime I corresponds to two branches). The point 0 bn corresponds to regime III (or branch 3) and 
the weights are 


Pen = 

= 

"Ten = 


0.275899, 


2 COS — 

V16 


1 H-= tan 

V2 


(e) 


0.446933, 


1 (^2 + 72 +^ 


2 + V2 


2.630986. 


(9) 


The other three solutions can be shown to correspond to a point in the dense phase (above the blue 
surface in figure]^; a point at the Ising ordering transition (on the blue surface); and a point in the dilute 
phase (on the red surface in figure 

The point 0ds is derived from another—simpler—integrable model in which only pg and pg are non-zero. 
We denote the loop weight in this completely packed model by n. It is equivalent to the Q = state critical 
Potts model [^, and is referred to as branch 0 in 17 . One can shown that it is equivalent to a model of 
the type 0 with non-zero pi via a simple trick |17| that we now recall. Indeed, consider a loop model with 


only vertices pg and pg and loop weight n = n+1. Draw each loop independently in black and white colour, 
giving weight n to black loops and 1 to white loops. The trivial weight given to white loops means that we 
do not need to keep track of their number, so the colour white can be considered invisible. This produces a 
model of the type Q in which pe = P 7 = 0. Starting from the completely packed model with normalised, 
isotropic weights, pg = pg = 1, we obtain a dilute model with pg = Ps = P 4 = Ps = 1 and pi = 2. Setting in 
particular n = 0 we finally obtain 


PDS = 0, 




"rDS = 2 . 


( 10 ) 


By means of the h = n -I- 1 mapping, this model is equivalent to percolation (the Q = 1 state Potts model), 
just as the theta point model originally exhibited in 12 . In particular, the critical exponents at 0ds coincide 


with those of Q. 


3 Relation with the 0(n) loop and Izergin-Korepin vertex models 

An important fact—which is at the root of the exact solution in [2Q| — is that the point 0bn is just the 
particular n = 0 case of a more general solvable loop model related to the Izergin-Korepin ^-matrix. We 
now review this in some detail. 

®For completeness, we notice that there is another point in the dilute phase which is extremely well studied and whose 
position is conjecturally exactly known. It is the critical point Kc of the combinatorial SAW model, with p = 1 and r = 0, 
in which vertex-avoidi ng w alks are simply weighted by K per monomer. This is well suited for exact enu mera tion studies, of 
which the most recent [38] confirms to almost 12 decimal places the validity of a twenty year old conjecture |39| that Kc is the 
positive real root of the polynomial 581/^"^ + 7K^ — 13. 
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(‘2\ 

The Izergin-Korepin i?-matrix (also known as the ' i?-matrix) is a solution of the Yang-Baxter equa¬ 
tions acting on two copies of the spin-1 representation of Uq{sl 2 )- It is most conveniently expressed in terms 
of the quantum group projectors Pg on spins s = 0,1, 2: 


oc Pa + 




1 


— X 




( 11 ) 


40 


We will set in the following q = e'^ and x = in order to match the notations of Ref. 

While the most natural context to interpret such an P-matrix is in terms of a 19-vertex model, it turns 
out that a few simple transformations allow a reinterpretation in terms of a dilute loop model 41 ^ 

In order to discuss this—and in view of the large literature on the subject—it is useful to review nor¬ 
malisations and conve ntio ns in detail. First, le t us fix the normalisation of the P-matrix. We start from the 
matrix Pgm given in [^, which is related to (|ll[) by 


_ (q^-x) {q^+x) 

-^GM ^ 5 '' 5 


xg2 


where q = We then define our model by 


P = 


(sin (i) - sin (1^)) 


Pgm ■ 


( 12 ) 


(13) 


Next, we graphically represent the edges with S'z = 0 by empty edges, and those with Sz = ±1 by 
oriented lines. We then apply the following gauge transformations on the matrix elements of P: 



with a = ie'^ Then perform the following transformation 



with /3 = —ie Considering now the model as a gas of dilute loops, each loop carrying the sum over its 
two possible orientations, the local configurations around a vertex are precisely those shown in 0. The 
latter gauge transformation amounts to assigning to every closed loop a weight 


n = = —2 cos( 27 ) , 

and the weights pi,..., pg assigned to each vertex configuration are the following 

sinuzB sin(3AzB — uzb) 


(14) 


Pi = 1 + ■ 


P2 = P3 = 


Pi = P5 = 


P6 — Pi — 


sin 2 Azb sin 3 Azb 
sin(3AzB — Uzb) 
sin3AzB 
sin MzB 


sin3AzB 

sin Uzb sin(3AzB — uzb) 


Ps = 


P9 = 


sin 2Azb sin 3Azb 
sin(2AzB — Uzb) sin(3AzB — uzb) 
sin 2Azb sin 3Azb 
sinuzB sin(AzB — uzb) 
sin 2Azb sin 3Azb ' 

—2 cos 4Azb • 


(15) 


'This is somewhat similar to what happens for the 
with a completely packed 0(n) loop model [42]. 


( 1 ) 


model .R-matrix, which is both related with the 6-vertex model and 
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The model Q with weights (151 is precisely the integrable dilute 0(n) model [T7}|^ , here written in the 
conventions of In the following we parameterise the anisotropy by A (spectral parameter) and the loop 

weight by 7 (crossing parameter), as in (14). These are linked to the notation through the relations 


TT 7 

Azb = 2-2 


UZB 


= iA. 


(16) 


In all the following we also rescale p 4 and ps by a factor —1. As noticed in 17 this simple gauge 


transformation changes neither the partition function nor the correlation functions, and it has the advantage 
of rendering all the Boltzmann weights at the polymer point positive, as shown explicitly in 0 - This is an 
important point to us, since we want to show that this dilute polymer model has a non-compact continuum 
limit (i.e., with a continuous spectrum of critical exponents), although the only non-unitary ingredient stems 
from the non-locality of the loop weight n. This should be compared in particular with the earlier examples 
of statistical mechanics models with a non-compact continuum limit, that were defined in terms of negative 


Boltzmann weights 26 27 29 


With these conventions, the weights (15) admit two isotropic points which are 




(17) 


and for each of those there are two polymer (n = 0 ) points in the fu nda mental interval 7 € [ 0 , tt ], na mely 
7 = f and 7 = ^ ■ The correspondence with the three regimes of Ref. 


20 


and the four branches of Ref. 


17 


is as follows. With u = m +, we are in regime I, corresponding to dense polymers (branch 1) for 7 = 7 , anc 
to dilute polymers (branch 2) for 7 = Explicitly this leads to the weights 


Pden — 
^den — 

'^den — 


y2cos — « 1.38704, 
16 


y/2 cos ^ — 2 sin ^ 



1.00315, 


0.783227 


(18) 


and 


Pdii = 


ATdii = 


T'dil = 


/- . 37r 

V 2 sin — 
16 

y2sinf| 
1 -I- cos I 


0.785695. 


0.408391, 



0.675577. 


(19) 


With u = U-, we are in regime III (branch 3) for 7 = 7 , corresponding to the point 0bn already given in 
91), and in regime II (branch 4) for 7 = ^ with weights 


Pii 


r- 37r 

= v 2 cos — 
16 


1.17588, 


Kn = 


TII 


4 sin 77 sin 7 


1 - 


) sin I sin^ 77 


15.4476, 



-0.0897902. 


( 20 ) 


Unlike the other integrable points 0 , ( |l^ and ( |l^ , this latter point ( [^ does not seem to allow for a gauge 
transformation that will render all Boltzmann weights explicitly positive. 


^The subscript ZB refers to Zhou and Batchelor 4l]. 
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4 Transfer matrix spectrum and critical exponents 

4.1 0(n) model transfer matrix 

Consider the 0(n) loop model on a square lattice of size Lx H sites in the horizontal and vertical directions 
respectively, and take periodic boundary conditiorrs in the former. The corresponding (grand canonical) 
partition function can be written as 


Zl.h 


c 


( 21 ) 


where the sum is over all possible loop configurations, Ni (with i = 1,..., 9) is the number of occurencies 
of each type of vertex Q, and Moops is number of closed loops within each configuration. Note that both 
closed contractible loops and loops winding around the horizontal periodic boundary condition contribute to 
Moops- It is possible to give a different weight to these two loop types by twisting the model, and we shall 
make use of this possibility in the sequel. 

It is convenient to introduce a transfer matrix formalism 


17 in which (21) is computed by decomposing 


the lattice into horizontal slices. Each configuration of the vertices within the lowermost t horizontal slices 
induces a connectivity among the L edge mid-points that form the intersection of the lattice with a horizontal 
line cutting the system between row t and row t + 1. Some of these L points can be empty, and the remaining 
points are either connected pairwise via an arc (a contiguous part of loop situated below the intersection 
line), or they connect all the way through the system via a leg (a contiguous loop segment that sends at the 
bottom of the lattice). 

The set of all possible connectivities for a system of size L = 3, labeled by the number of legs £, can be 
drawn as follows: 


o 

II 

£ 

1 

= 1 

- II 

. to 

£ = 3 

1 1 1 

V/ • 


1 ■ 

- 1 1 

1 1 1 


1 

• 1 

1 ■ 1 




1 \ 




^ I 


We note that since the vertices ([^ do not allow the crossing of loop segments, the configurations can all be 
drawn in a non-crossing fashion. In particular a leg is never trapped inside an arc. The horizontal periodic 
boundary condition however implies that for £ > 0 some of the arcs will need to cross the boundary. 

From a configuration corresponding to a given set of connectivities at height i7, the possible configurations 
at height H + 1 are obtained by acting with all the compatible vertices ([^ on each of the L vertices within a 
row. Each connectivity carries a statistical weight, corresponding to the Boltzmann weights appearing in the 
partition function (21) for a partially constructed lattice. The addition of the {H + l)th horizontal row can 
thus be represented by the so-called geometrical (row-to-row) transfer matrix , whose non-zero elements 
provide the transitions between configurations at height H and H + 1. More precisely, the matrix element 
is equal to the product of Boltzmann weights (j^ for the L vertices within the {H -|- l)th row—including 
a factor of n for each closed loop—summed over the diagrams from 0 that produce connectivity j for the 
system of size H + 1, starting from the connectivity i for the system of size H. The transfer matrix can be 
depicted as follows 


1 2 


L 











where each crossing represents the i?-matrix of section]^ in its loop representation, and it is understood that 
the left and right ends of the horizontal line have to be joined to recover the periodic boundary conditions. 
This is exactly equivalent to equation (2.12) in Ref. [l7] . 

It is easily seen that the transfer matrix has a block-triangular structure, as the number of legs i 
can only be lowered (by contracting one leg with another) or kept fixed from one row to the next. As will 
become clear soon all the information we need is encoded in the transfer matrix’s eigenspectrum, so we can 
focus on the diagonal part, and hence study separately the different subsectors of fixed i. Note that the 
transfer matrix is not symmetric, which means it is not necessarily diagonalisable. It could also be that some 
of the eigenvalues are complex. However, numerical study shows that, for generic n, the matrix is indeed 
diagonalisable and all the low-lying eigenvalues (which are those that matter in the scaling limit) are real. For 
degenerate values of n (such that g is a root of unity), the transfer matrix is usually not fully diagonalisable, 
and exhibits Jordan cells [43[|44| , even in the low lying sectors, indicating that the corresponding conformal 
field theory is logarithmic [45j7This is, however, a technical detail which does not matter much here, and in 
general we do not distinguish genuine and generalised eigenvalues. 

Since the model is critical, this geometrical transfer matrix has a large number of very close eigenvalues 
Aa , whose scaling with the size of the system is related with critical exponents, following the usual conformal 
invariance prediction. Calling the largest of these, and A^^^ any of the following eigenvaluesj^ we 


have 46 47 


log A 


(i) 


logAl,^) 


logA[,^^ 


f , fr- 

/oo-^+o(L 


2ttXc 

~1X 


+ o{L ^). 


( 22 ) 


where fao is the free energy per vertex of the infinite system, c is called the central charge, and Xa = Aq,-|-Aq, 
where (Aq, Aq.) are the holomorphic and antiholomorphic conformal weights of the operator associated with 
the state a. The central charge and set of conformal weights associated with the leading transfer matrix 
eigenlevels encode the whole set of critical exponents, as will be illustrated in the following sections. 


4.2 Relation with the IK transfer matrix, and exact results 

In order to investigate the critical content at the ©bn point, our next task is therefore to identify the leading 
eigenlevels of the 0(n) model transfer matrix, or alternatively, thos e o f the Izergin-Korepin model transfer 
matrix. The latter is defined similarly in terms of the matrix (II), and it acts in the product (C^) 
of spin-I representations (see 28 for more details). It commutes witn the total magnetisation 




(23) 


which is the vertex model counterpart of the geometrical quantum number £ mentioned earlier (see also 
below). A crucial point is to control the relationship between the loop model and the different gauge 
transformed variants of the vertex (Izergin-Korepin) model. There are well-defined procedures to do so, based 
either on algebraic considerations [^, or on explicit calculations of generating functions of levels 51 . Note 


that this is not only a technical question: as the loop model is essentially non-local, the kind of question— 
typically, a geometrical correlation—one wants to answer will affect the details, or even the nature, of the 
correspondence. In what follows, we restrict ourselves to the eigenvalues of the ‘natural’ polymer transfer 
matrix, discussed in 17 : it encodes in particular all the exponents we are interested in, and that have been 


studied in the literature. As a matter of fact, the precise relationship between the 0(n) and Izergin-Korepin 
eigenvalues depends crucially depends on the number of through-lines £: 


®It is convenient to view our two-dimensional statistical mechanics model as a one-dimensional quantum model (spin chain) 
evolving in imaginary time. In this picture the leading and next-to-leading transfer matrix eigenlevels are the ground state and 
low-lying excitations of the corresponding quantum Hamiltonian. We shall frequently refer to this equivalent point of view in 
the following. 
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• In sectors £ > 0, the eigenvalues of the loop transfer matrix are found in the spectrum of the vertex 
transfer matrix with periodic boundary conditions, 


= Tr^ [Rhi^h2 ■ ■ ■ 


hL 


(24) 


in the sector of magnetisation m = i (or, equivalently, m = 
horizontal space h = 


= —£). Here the trace is taken over the 


• Things are slightly different in the sector i = 0 where the loop model allows for noncontractible loops, 
that is, closed loops that wrap around the cylinder. These must be weighted n —the same as the 
contractible loops—which requires taking for the vertex model the following twisted transfer matrix 




(L) 


vertex 


Tr,,( 


d( 2) p(2) 

^hl ^h2 




(25) 


in the sector of zero magnetisation, m = 0. This transfer matrix differs from the periodic one by the 
presence of a boundary twist term , where is now the magnetisation along the horizontal (or 
auxiliary) space h. Summing over the two possible orientations of a non-contractible loop yields a 
weight e*'^ + e”"^, which can be made equal to n of Eq. (14) by taking Lp = tt — 27 . The eigenvalues of 


the loop model transfer matrix are a subset of those of the vertex model with this particular value of 
the twist. 


Let us again recall that the Izergin-Korepin model exhibits several regimes. We find it convenient, in 
order to identify these regimes, to restrict to 7 G [ 0 , 7 r], and use the symmetries of the weights. The model 
( 11 ) has two particular values of the spectral parameter, x = ii?®, for which the weights are isotropic. 


The case x = iq 2 corresponds to the so-called regime I, whereas the case x = —iq^ splits into regime II for 
I < 7 < TT, and regime III for 0 < 7 < The point Qbn that we are interested in corresponds to a; = —iqs 


and 7 = f • In the notations of [^, this is'k=|,0 = —Regime III has been studied in details in 
and we will heavily borrow results from this reference in the following. 


28 


4.2.1 Structure of the leading transfer matrix eigenlevels 

Regardless of whether we consider the vertex or loop formulation, the transfer matrix commutes with the 
lattice momentum operator which is defined so that e'^^ amounts to a horizontal translation by 

one lattice unit. The transfer matrix eigenvalues can therefore be classified according to the value of the 
momentum, with p = 0,1,..., L — 1, which in turn is equal to the conformal spin, A — A = 

From the integrability of the Izergin-Korepin model, the transfer matrix eigenvalues can be computed 
exactly in terms of a set of complex numbers, the so-called Bethe roots which are solutions of the Bethe 
Crucially, the number of Bethe roots (and equations) increases linearly with the size L 


Ansatz equations 28 


of the system (more precisely, the number of Bethe roots associated with states in the sector of magnetisation 
m is mi = L — to), whereas straight-forward diagonalisation of the transfer matrix would involve solving a 
set of equations whose number increases exponentially with L. 

It is convenient to discuss first the structure of the leading transfer matrix eigenlevels in regime III for 
the untwisted case (i.e., with periodic boundary conditions). The ground state (largest eigenvalue) is found 
in the to = 0 sector, and corresponds to a root configuration made of a ‘Fermi sea’ of j 2-strings (pairs of 
conjugate roots with imaginary parts close to±(j — j)). The leading excitations with zero momentum in 
the TO = 0 sector are obtained by taking an arbitrary number j of 2-strings off the Fermi sea and replacing 
them by j pairs of anticonjugate roots on the axis of imaginary part The story is essentially the same in 


^®Note that the loop model in the sector £ > 0 has more eigenvalues than the vertex model with periodic boundary conditions 
and magnetisation m = iL£. 

This is because, for a given £, one can take each of the £ through-lines around the axis of the cylinder without affecting the 
Boltzmann weights, if every such line acquires in doing so a phase which is an root of unity. This means that some of the 
polymer eigenvalues are obtained by taking the vertex model with m = £ and twisted boundary conditions, the twist being of 

2i7rp _ 

the form e ^ , with p A i = 1 and p integer. The corresponding critical exponents then satisfy A — A = p + integer. These 

exponents are crucial in determining for instance the winding angle distribution of the polymers [52] . They are, however, not 
so important for our purposes, and therefore we will not study them further here. 
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other magnetisation sectors, and from now on we will therefore label the leading eigenvalues of the transfer 
matrix by two integers, (m,j). 

We can now extend this to generic values of the twist. There are minor changes in the qualitative 
description of Bethe root configurations, but the global structure of eigenlevels unchanged. Going back to 
( 22 ) we therefore have Aq = A(o q), and we shall define the critical exponents Xmj = (where 

we recall that in the zero-momentum sector) or, equivalently, the effective central charges 

= C 12x772 = Cq^O 12 X 777 j'. 


4.2.2 Sectors with £ > 0 (m 7 ^ 0) 

In the case of non-zero magnetisation, we have already seen that the twist is zero. It follows from our Bethe 
Ansatz analysis that the leading effective central charges have the form 


12 


+ (JV„,,7 


Ain) 


Att 


[Bm,jil) + '^0g LY 


where j labels the excitations. Here, the Nm,j are integers satisfying 

1 


and the function A( 7 ) is 


= 2j + -(3-(-1)-) 




(26) 


(27) 


(28) 


The functions Bmji'j) are not known accurately at this stage, but are believed to be universal functions 
with no L dependence at this order. 

In the limit L —>• 00 , the L-dependent term on the right-hand side of (26) goes to zero at fixed j, suggesting 


that the ground state in this sector is infinitely degenerate. Although this is, in a sense, true, it is better 
to interpret what happens by observing that, if one scales j with logL so as to keep s ^ finite (and 


continuous), Eq. (26) can now be interpreted as the signature of a continuous spectrum of critical exponents. 


Such a spectrum is not so familiar in statistical mechanics. It occurs frequently in string theory, where the 
‘targets’ of associated conformal field theories are not compact. This non-compactness leads to continuous 
spectra of critical exponents just like, in ordinary quantum mechanics, free particles on a non-compact space 
(e.g., a real line) have a continuous spectrum of eigenvalues. In most problems of statistical mechanics, by 
contrast, these ‘targets’ are compact. This is well-known in the case of the Coulomb gas representation, for 
instance, where the bosonic field is compact, that is, d) = <& -|- 27ri?, where i? is known as the ‘compact- 
ification radius’ 53 . In the IK model the continuous spectrum can be interpreted as arising instead from 


a ‘non-compact boson’, that is, a bosonic degree of freedom $ for which, formally, i? = 00 . More precisely. 


the log L dependence in Eq. (26) can be thought of as arising from an effective compactification radius Ji(L) 
that diverges as L 


00 . 


We represent qualitatively the spectrum of critical exponents in the sectors with m > 0 as shown on the 
left panel of figure Each spectrum has a ‘socle’, above which a continuum starts immediately. We will 
discuss below what this means for the correlation functions. 


4.2.3 The sector with no through lines (to = 0) 

The main difference with the case of £, to 7 ^ 0 is that the leading eigenvalue (which determines the central 
charge) is separated from a continuum of critical exponents by a gap. This has to do with the fact that, 
in the IK model, the whole spectrum contains both continuous and discrete states —^just like, for quantum 
mechanics on the line in the presence of a potential, one can have a spectrum made of both a continuous 
and a discrete part. The ground state (with j = 0 and to = 0) at this value of the twist is given by a discrete 
state, and leads to the central charge 

1272 

Cq.O — —1 H--r • 

7r(7r — 7 ) 


(29) 
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X = A + A. 



X = A + A. 



^0 


J = 0 


m = 0 


Figure 3: Structure of the spectrum of critical exponents in the sectors m > 0 (left), and m = 0 (right). 
In the latter case there is a gap be twee n the gro und state (j = 0) and the continuum of j > 0 excitations, 
whose magnitude is obtained from (29) and (30) as Xg = — - 


_ C0,1 —co,o _ 7 

12 TT —7 


This coincides with the prediction of Nienhuis et al: see Eq. (7.3c) in 20 , where the notation 9 = —7 is 


used. The excited states j > 0 in the m = 0 sector form part of the continuum, leading to the effective 
central charges 

co^- = 2-^-12(Aogf , (30) 

' [Bo,, (7) +log 

These observations can be summarised as shown in the right panel of figure 

From all this we can deduce the so-called watermelon exponents of the loop model. The m-leg exponent 
is defined as 

Cm,0 Co,0 


— 


12 


(31) 


where the effective central charge is measured at the appropriate twist parameter as detailed above, and the 
L —> 00 limit is taken such that the non-compact part does not contribute. We find then 


1 m ^7 


7 


7 r( 7 r — 7 ) ’ 


(32) 


which agrees, once again, with 20 : see Eq. (7.5) in that reference. 


4.2.4 Leading exponents at the polymer point 

From these results we recover the central charge c = 0 and the leading watermelon exponents at the polymer 
point, which corresponds to 7 = j, hence loop weight n = 0 : 


— 


m 

16 


(33) 


This means in particular that 


^^=-48’ 


^^ = I 2 


^4 = X ■ 
0 


(34) 


This leads, by the usual scaling relations |3| (namely ry = 2xi, z/ = 2 -X 2 ’ ^ = 2 — ry) to the exponents 

given in (§. In addition we identify, like M the usual theta point, the four-leg operator with the thermal 
operator coupled to the monomer-monomer attraction energy, which drives the problem away from the 
(tri)critical point towards the dilute or massive phase. By the scaling relation X 4 = 2 — ^ one finds z/' = |, 
and thus the crossover exponent 
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5 Physical consequences of non-compactness at the 0 bn point 


5.1 Correlation functions 

In the polymer problem, one is typically interested in the correlation functions of observables which are most 
easily defined in lattice terms, such as the polymer limit (n —^ 0) of the spin-spin correlation function. The 
latter can be expressed as a sum over all configurations of (in this case, vertex-interacting) self-avoiding 
walks joining two given points: 


Gi(J, J) = 


E 

VISAW: /->,/ 


K’ 


umber of monomers number of doubly visited site: 


umber of straight segments 


(36) 


Such correlation function become, at the critical point and in the continuum limit—that is, here, at distances 
much larger than the lattice spacing—a complicated sum over correlation functions of conformal fields. The 
point is that a given lattice observable never corresponds exactly to a pure scaling field: rather, it can be 
expanded as an infinite sum over scaling helds, typically of the form 




(37) 


where the Oi are pure scaling fields whose two-point function is normalised to l/{zz)^' (for simplicity, we 
only mention scalar contributions in this analysis). The analysis of the correlation function (361 must take 
into account both this expansion of operators and a similar expansion of the lattice Hamiltonian, which 
differs from its continuum limit counterpart also by a sum of irrelevant operators. In models with a discrete 
spectrum of critical exponents, this simply leads to 


Gi(J,J)cx ( 


G(^) 

JiJ 


2x'-, 


(38) 


where x'l > xi is the exponent of the first corrections to scaling term. The discrete spectrum here guarantees 
that the contributions of the various correction terms to Gi can be well separated at large distance, and that 
the asymptotic behaviour is fully determined, at leading order, by xi. This problem is well-known, and has 


been studied extensively, for instance for the six-vertex model in 54 


Meanwhile, when the theory admits a continuous spectrum of critical exponents, one expects sums such 


as (37)-(38) to be replaced by integrals: 


Gi(/,J) 


ds 


fljs) 

2xi+2ais^ 


We have here written the continuous part of the spectrum simply as xi -\- ois^, with s 


(39) 


l4x and ai = 


4H(7r/4), by expanding (26|-(27) to leading order. We could absorb ai into the definition of the quantum 


number s, but we refrain from doing so, because s has a precise meaning when the conformal field theory is 
interpreted in terms of a non-compact coset sigma model. The amplitude /i(s) is meanwhile determined by 
several factors, including the matrix elements between the lattice spin observable and the eigenstates with 
quantum numbers (m,j) = (l,z) and the (properly defined) density of states. Even though the theory is 
non-unitary, we expect the function fi to be positive^] 

Interestingly, expressions such as ( [M| have appeared before in statistical mechanics. Most noticeably, 
at the plateau transition in the quantum Hall effect for instance, it is expected that the g**' moment of the 
transmission coefficient between two point contacts at distance r reads 55 


(n= J \{Vq,V*0\Xq)\\-^^^fi{X)dX 


where the scalar products and the measure /i(A) are known, while the exponents only have a conjectured 
form. In this case as well, the presence of the integral is ultimately related with the non-compactness of the 
target. 


'^^This positivity was checked numerically on cylinders of sizes L = 4,6,8,10, where we explicitly computed the matrix 
elements (l,i| Z;f=i 0>)- 
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The practical point of an expression such as (39) is that the term going as r —which is the only 


one identified by the naive determination of the one-leg operator—predicts the leading behaviour of the 
correlation function only at extremely large r. It is difficult to say much more without having a better idea 
of the function /i, but it is nonetheless clear from the outset, that any attempt to match an expression such 


as (39) with a pure power law r will lead to an effective exponent x > Xi. In the specific context of 


numerical simulations, it is of course well-known that the results based on the study of a system of linear size 
L provide only an effective exponent xi{L), which will have to be extrapolated by using finite-size scaling. 
Our point is that such extrapolations are usually thought to contain only power law corrections, and this 


does not correctly take into account the presence of logL in (26). 
Writing now 

poo 

/o 


Gi(r) =r 


— 2xi 


(40) 


there are known techniques 56 to extract an asymptotic expansion for Gi based on the analyticity properties 


of /i near the origin. We will restrict ourselves to the leading term here, and simply observe that if 


/i(s) — s* [6o + s^i + • ■ •] 


for 


0 , 


(41) 


we will have 


Gi(r) = r-2^i(lnr)-(^+^)/2 


Co 


Cl 


-k.. 


(42) 


5.2 Grand canonical Monte-Carlo simulations 


Let us compare these predictions with Monte-Carlo sampling of the function Gi(r), both at the point ©bn 
and at the point ©ds- the latter point the peculiar effects of non-compactness are not expected, and 
indeed it is well-known 12 that X\ = ^. The polymer configurations are sampled on square lattices of sizes 


L X L with periodic boundary conditions in both directions, with one polymer end fixed at the center. The 
second end evolves according to an improved version of the backbite algorithm 
into account the Boltzmann weight of doubly visited sites and straight segments 


57 58 that properly takes 

0 

The trivial, empty configuration which is just represented as a point-like polymer at the center of the 
L X L lattice is assigned the winding numbers Wx = 0, Wy = 0. Whenever the polymer crosses one vertical 
or horizontal boundary the corresponding winding number is increased by ±1. We proceed as follows: 


1. Start from the trivial, point-like configuration; 

2. Run 500 x L x L time steps to decorrelate from this initial configuration; 

3. Then run 10^ x L x L time steps, and sample the configurations with Wx = Wy = 0, measuring the 
end-to-end distance r between the two ends of the polymer. 

This produces a histogram of of the form 


H{r)= Y. Gi(0,J) + ..., 

|J-0|2=r2 


(43) 


where the dots indicate the finite-size corrections related to the fact that the polymer interacts with itself 
if part of it goes around the torus in either direction. In practice, we produce for each size several such 
histograms (between 10 or 20 depending on the size), and compute the corresponding average and error 
bars. The function Gi(0, J) can be extracted from H{r), either by a binning procedure, or for each integer 
value of by dividing H(r) by the number of different possible J such that | J — Op = r^, eliminating those 
for which no such J exist (for instance = 3 cannot be obtained on the square lattice). 

be precise, the pure backbite algorithm [57| applies to the fully-packed (Hamiltonian walk) case. It can be adapted 
to dilute walks by adjoining grow and retract moves, in a Metropolis scheme, following the last few lines of [58] . Finally, the 
possibility of having doubly visited sites introduces a few additional complications. The most subtle of these concerns the 
situation where a polymer end point resides precisely at a doubly visited site, in which case a special backbite-type move is 
introduced, under which the end point moves infinitesimally across the other loop strand visiting the concerned site |59|. 
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Log L 


Log G(L/4) 


- 0.5 - 

Log G = - 2 Xi Log L + A - B / L 



Log G(L/4) 



Figure 4: Fit of logGi (0.25L) against logL, measured from Monte-Carlo sampling at the points ©ds (top) 
and ©BN (bottom) for L = 10,..., 100 (the error bars are also plotted, but are smaller than the size of the 
symbols). In the first case, a fit of the form (441 seems in good agreement with the numerical data, and 
yields xi ~ 0.245. In the second case, such a fit is clearly not appropriate (orange curve), and we find much 
better agreement with (451 (blue curve). 


After normalizing Gi(0,1) = 1, for a separation of one lattice spacing, the function we eventually consider 
is Gi(r) = Gi(0,r), measured for r equal to a fixed ratio of L, that is, r = aL with a ^ 1. In practice 
a = 0.25 will turn out convenient. This ensured that with respect to the total size of the lattice the polymer 
is almost closing on itself, and therefore the corrections coming from walks going around the torus and back 
should scale as where X 2 is the two-leg watermelon exponent, which is positive for both ©bn and 

©BS. From the previous section we expect for large L: 

• For Qds (and other usual cases) 

logGi(aL) =-2a;i logL -I- A -I- G ^ -I- ... (44) 

• For ©BAT, 

log Gi (aL) = -2xilogL + A - log log L + + ■. ■ (45) 

2 i/log L 

As shown in figure (where in the former case we set G t , which takes correct 

account of the corrections to leading order, while in the latter case B = '^^), these expectations are nicely 
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corroborated by our results. At 0 ds we find xi ~ 0.245, which is close to the expected xi = j. At 0 bn 


we used a fit keeping only the first three terms on the right-hand side of (451, and we plugged in the known 


value of xi = — This yields z ~ 1.7; but note that this estimation is very sensitive to the presence of 
higher corrections in the fit, and therefore we shall not trust it too much quantitatively. It is quite striking 
that, even though Xi < 0—and therefore at very large distances the function Gi(r) should be an increasing 
function of r —, at the distances that could be achieved in these simulations Gi(r) is a decreasing function 
of r, leading to an effective exponent cci > 0 at intermediate distances. 

Albeit the numerical results presented in this section seem rather satisfactory, it could of course well be 


that /i has a behaviour quite different from that predicted by (411. Unfortunately, at this stage very little 
is known about this problem. 


5.3 Further speculations 


In most of the lattice Monte Carlo simulations presented in the literature 22 
is however the canonical (fixed length) partition function, 


the quantity usually studied 




umber of doubly visited site: 


number of straight segn 


(46) 


VISAW: I^J 


Even if the function fi were known, going from Gi (I, J) to Zi (/, J) would require making hypotheses on 
the scaling form of the correlation function in the vicinity of the critical point, K ^ Kbn, and inverting 
the corresponding expressions of Gi considered as a Laplace transform of Zi. There seems little point in 
speculating about what may happen until more detail is available for the function /i. A short-cut is however 
provided by the following argument. Recall the scaling relation 


7 


= 2 - 77 , 


(47) 


where 77 = 2xi. Since the two-point function involves a mix with values of Xi greater than the bottom of 
the tower, this should lead to estimates {v being what it is) for 7 which are smaller than the true value. 
Indeed, the numerical result 22 reads = 1.045 < 1.152. Meanwhile, from a ;2 = 2 — ^ we see that 


2 — X2 ' 


Since similarly X 2 is decorated by the tower on top of the minimum value, it will appear larger 
than its asymptotic value as the polymer length goes to infinity, which means that 2 — a ;2 will appear smaller 
and thus ly larger. Indeed, the numerical result 22 is = 0.576 > 0.522. So the discrepancy between 


numerics and exact results is compatible with our scenario. 


6 The polymer phase diagram 

The main question now concerns the relationship between the theta points 0bn and 0dSj and more generally 
the nature of the phase diagram for the VISAW. 

It is certainly tempting to consider the 0bn point as an unphysical, infinitely multi-critical point. This 
is because the spectrum in the sector with no through-lines is continuous above the main gap given by 

This means there is an infinity of relevant operators in the thermal sector, with 
A < 1, corresponding formally to an infinity of relevant directions, hence an infinite order of criticality. 
On the other hand, it is troubling that the the 0bn point seems so natural, and only necessitates the fine- 
tuning of a couple of parameters in our phase diagram. Moreover, our experience with the stability of RG 
fixed point with continuous spectra of relevant exponents is limited: one should not be too hasty in deciding 
that the 0bn universality class is unphysical. We must now try to understand this universality class and its 
vicinity in more detail. 


6.1 The second virial coefficient 


At 0BN (7 = f )i the ground state of the Izergin-Korepin model with twist <75 is determined by two possible 
values of the central charge: 


c((p) 


c* = 2 - , for (/? < f , 

-I + , for(/ 7 >f. 


(48) 
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We must take v? = | to make the weight n = 2 cos ip of non-contractible loops equal to zero, which leads to 
c = 0. This is indeed the value expected for the theta point. It is now interesting to consider what happens 
near ip = ^. First, the leading eigenvalue of the transfer matrix takes the form 


logA^^^ 

L 


foo 


TTC^p) 

6L2 


+ o{L , 


(49) 


where fao is independent of p, and there are no corrections of order A because we use periodic boundary 
conditions. When p ^ loops that wind around the axis of the cylinder are allowed, with fugacity 
n = 2cos(^. Setting ‘/J = f ~ f > this becomes fi = e + 0{e^). Meanwhile, we can expand the central charge 
at this order, and thus the logarithm of the leading eigenvalue: 


logA[,^) = 


1 

L 



(50) 


where we have used the natural normalisation where the partition function at n = n = 0 is unity (hence 
/oo = 0). Imagine now taking a long cylinder of length L': the partition function to leading order will be 



Meanwhile, this partition function can be expanded in powers of n: 


InZ = fizi + 



(51) 


(52) 


where zi is the partition function with one non-contractible loop anywhere on the cylinder, and z^ the one 
with two such loops. Comparing, we find the result 



J_ 
L Gtt 


(53) 


The left hand side combination can be considered as a sort of second virial coefficient . Such terms are not 
particularly inspiring in general. What is remarkable here is that this difference is positive. This is related 
to the positive curvature of the second line of (48) for the central charge as a function of p. In models of 
polymers, such a term is usually negative, like in the first line of (48), leading to the opposite sign for the 
combination of partition functions (53). It is easy to understand why such a term is negative for instance 
for ordinary self avoiding walks: the point is that every term in Z 2 occurs also in ^zf, while there are plenty 
of terms in that are not in Z 2 —all the terms where two non-contractible loop from each of the zi have 
a non-zero overlap. When the model involves bending and attraction, it is not so easy to establish the sign 
of this term before hand. For instance, there are terms in Z 2 where two non-contractible loops are close to 
each other, contributing a mutual attraction energy which is not present in the corresponding (geometrically 
identical) term subtracted in ^zf. Nonetheless, in all models known to us—in particular the theta point 
©DS—the combination remains negative. The fact that it is positive in our case indicates the attraction 
between monomers at the ©bn point is unusually strong, a significant fact confirming the belief that ©bn is 
a higher order critical point. 


6.2 Relationship with the dense phase 

To proceed, we make the simple but crucial observation that the radius of the compact boson for a given 
value of 7 —that is, for a given value of the fugacity of contractible loops—in regime III is the same as the 
radius for the compact boson describing the corresponding dense phase of the vertex or 0(n) model. The 
dense phase of the vertex model is obtained by increasing the fugacity K of the monomers in the loop version 
of the problem, and is in the same universality class as the “completely packed” model, which is nothing 
but the ordinary six-vertex model. The corresponding dense phase of the 0(n) model is obtained in the 
same way, the only difference being that non-contractible loops get fugacity n, instead of two for the vertex 
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model. This observation can be illustrated by considering the polymer case (n = 0). The m-leg (watermelon) 


exponents are found to be 28 

while they are 


1 

(54) 

■ It ~ 6 ’ 

— 4 

(55) 

16 

are quite close, since 



(56) 


The identical term arises from the fact that both theories involve the same compact boson. The shift of 
— ^ occurs because of the non-compact boson. More precisely, the scaling of the associated eigenvalue in 
the dense polymer case 


log A, 


(i) 




L 


(57) 


involves 


^ - 


c 

12 


m 

It 


1 

12 


_ = - 




12 


with (P = —2 for dense polymers, while we have for our critical point 


- .-V — 


C 

12 


m 

Te 


- - = - 
6 ^ 


r.2B 


12 


Here, = 1 denotes the central charge of one boson, while = 2 is the central charge of a system 
of one compact and one non-compact boson—or the effective central charge of the Black Hole theory. We 
see therefore that if we took the leading eigenvalue in the m-leg sector and made the non-compact boson 
massive, the scaling would simply evolve from our model to the dense polymer case: 


log A, 


(i) 


L 


-foo+^^(x ^ 2 ) 


f + — 

Joo -I- ^2 




.IB 


12 


(58) 


Meanwhile, we also notice that the central charge (that is, the scaling for the ground state energy) is given 
in our model by a discrete state, which involves the non-compact boson. Recall indeed the general formula 


(48) in our case for a twist (p. We must take ^ to make the weight of non-contractible loops equal to 
zero, which leads to c = 0, while c* = —1. This result for c* comes from the central charge c = 1 of the 
non-compact boson, and the central charge of the twisted compact boson which must be c = —2. We see 
that, if we made the non-compact boson massive while leaving the compact one untouched, the discrete level 
would disappear, and the value of the central charge would become c = —2 as required for dense polymers. 
Therefore, in many respects, the point Obn — s,nd more generally, the critical point in regime III—looks like 
a dense polymer supplemented by a non-compact degree of freedom. 0 This confirms the peculiar nature of 
the ©BN critical point, whose features are profoundly different from those expected at an ordinary tricritical 
point. 


6.3 Probing degrees of freedom 

It is interesting now to study numerically the behaviour of the model when the monomer fugacity K is 
increased, K > ATbn = ©bn- In figurewe have plotted the average density of monomers per lattice edge, 
evaluated from the derivative of the (ground state) free energy /o with respect to the fugacity K (the results 
are presented both in the untwisted case where non-contractible loops are assigned a weight h = 2 and in the 
twisted case where all loops are given the same weight h = n). We find very clear evidence of a first order 

It should nevertheless be borne in mind that the two bosonic degrees of freedom are eventually coupled, leading to more 
subtle effects. 


18 










- L=4 

- L=6 

- L=8 

- L=10 

- L=12 


- L=4 

- L=6 

- L=8 

- L=10 


Figure 5: Average density of monomers — evaluated numerically from the ground state free energy 
/o = — j log Aq, as a function of the monomer fugacity K around the point 0bn- The top panel corresponds 
to the “loop case”, where both contractible and non-contractible loops are assigned the same weight n = 0. 
The bottom panel corresponds to the “untwisted case”, where non-contractible loops are assigned a weight 
n = 2. In the former case there is a level crossing at iF ~ ATbnj which is not present in the latter, but in 
both cases it is apparent that in the limit L —)■ oo the density of momomers is discontinous across ATbn- 
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Figure 6: Eigenlevels fi = log in the sector £ = 0, measured for A = 8 as a function of the monomer 
fugacity K around the point ©bat, in the “untwisted case”, where non-contractible loops are assigned a 
weight n = 2. The black continuous lines are the levels (0,j) (with j = 0,1,2 from bottom to top). For 
j 7 ^ 0, these levels become massive away from K = ATbn- Similar results are observed in other magnetisation 
sectors and for other values of the weight h. 


phase transition, with a discontinuity in the monomer density, therefore recovering the conclusions drawn 
from a corner transfer matrix analysis. In the K > IFbn phase, the fluctuations of density are not 


25 


critical, and the physics of the model is simply that of the six-vertex model, or ordinary dense loops. 

More generally, we refer to figure where it is apparent that the levels with j ^ 0 become massive on 
both the K < Kbn and K > Kbn sides of the transition point ©bn- 

These observations strongly suggest that the non-compact degree of freedom is associated with fluctu¬ 
ations of monomer density, since these fluctuations become non-critical precisely when the non-compact 
degree of freedom disappears. This is, after all, not so surprising. The emergence of a compact boson in 
the continuum limit of the six-vertex model follows from the reformulation in terms of a model of heights 
defined on the dual lattice, with heights on neighbouring sites defined recursively by the orientation of the 
arrows separating them. Study of the RG flow shows that the dynamics of these heights is described by a 
free bosonic theory. This theory is compactified because, on a torus, the periodicity of the heights cannot 
be guaranteed, and the boson much be identified with its shifts by 2ttR 53 . In dilute loop models, the 


corresponding vertex model (viz., the IK model) also admits edges that do not carry an arrow. In the regime 
I of the IK model, these edges do not have much implication on the continuum limit, as they simply lead 
to a renormalisation of the effective coupling constant (radius) of the free boson; the physics is essentially 
the same as the one of the dense model. It is however perfectly possible to imagine a scenario where the 
continuum limit of the theory is more complicated, and, while the arrow-carrying edges still correspond to 
height jumps in a height model with a dynamics described by a free boson at large scales, the arrow-free 
edges correspond to some additional non-trivial degrees of freedom in the continuum limit. This is in fact 
exactly what happens in regime II, where these edges are related with the emergence of an Ising degree of 
freedom in the continuum limit. 

What we propose now is that the edges that do not carry an arrow—and which correspond, in the loop 
model, to edges not carrying a monomer—have a dynamics described in the continuum limit by a non¬ 
compact boson. Note that, since the total number of edges on the lattice is of course a constant, the number 
of empty edges and number of monomer-carrying edges are related. This means that, in a certain sense, the 
dynamics of the density at the critical point of regime III should be related with the non-compact degree of 
freedom. 
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6.4 A model for K ^ K] 


BN 


It is now time to recall a few facts about the Black Hole sigma model which, according to our earlier work 28 


describes the critical 0{n) model in regime III. As discussed in the foregoing sections, the question of the 
twist plays an important role. The Black Hole sigma model as discussed, e.g., in 30 32,60 corresponds. 


strictly speaking, to the periodic vertex model underlying the 0(n) model. In the geometrical formulation in 
particular, this means that the non-contractible loops have weight two. The model where non-contractible 
loops have weight different from two is described by a twisted version of the sigma model. For simplicity, 
we start with the untwisted model, where we recall we have seen the same behaviour—in particular, the 
possibility of a first order phase transition—as in the twisted one. In the semi-classical limit, the action 
reads (for simplicity we do not write the dilaton term) 


47r 




1 + m 


(59) 


and describes a sigma model on a target with the shape of a cigar. The compact boson is obtained by 
considering the model in the limit of large |'I'| where the action can be approximated, setting = pe*®, by 


k 

47r 


d^x {d^9f 


(60) 


The non-compact degree of freedom is, very roughly speaking, the distance along the axis of the cigar, that is, 
Id^j. There are some crucial subtleties here having to do with normalisability, and the best way to understand 
them is to discuss briefly the mini-superspace approximation 60 . Physically, this is obtained by considering 


the model on a cylinder of very tiny radius, so that fluctuations in one direction can be neglected, and the 
model becomes effectively a zero-dimensional quantum mechanics problem. The Hamiltonian is then nothing 
but the Laplacian on the target, 




(coth r + tanh r) -\- coth^r d1\ 


(61) 


where we have set 'k = sinhre*®. There are no L^-normalisable eigenfunctions, but only ^-function normal- 
isable eigenfunctions. They depend on two parameters: one is n S Z, the angular momentum of rotations 
around the axis, and j = — ^ -I- fs is related with the momentum s along the r-direction (the axis) of the 
cigar. The corresponding eigenvalue of the Laplacian is 


a; = A -I- A = — 


2j'(j + 1) n 
k 


2 

2k 


(62) 


The full wave functions read (recall that in the mini-superspace approximation, the wave functions 
coincide with the expressions of the (conformal) fields in terms of the degrees of freedom in the action) 


r2(-j + M) 


And 


r(|n| + i)r(-2j-i) 

sinh'"' rF(j -k 1 -k 




(63) 


Setting j = — i -I- i I, and using the notation 
things in such a way that 


for the corresponding eigenfunction, it is possible normalise 


(CjC) =^nn' [27r(5(p-p')+ ^o(p',n)<5(p-kp')] 1 

where the scalar product uses the volume element dw = A: sinh 2r dr d0. 

One finds, for the classical reflection amplitude, the expression 


Ro{p,n) 


r(ap)r^ (i - f + §) 

(i + f + f) • 


(64) 


(65) 
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A crucial—and general—aspect of these wave functions is that they are not simple expressions of the 'h 
degree of freedom. Moreover, while these eigenfunctions have a ‘normal’ behaviour at small jd/j—namely 
(j>i^ oc for n > 0 (or for n < 0)—their behaviour at large 'b is definitely not usual in quantum field 
theory. To see this, it is convenient to decompose the eigenfunction as 


= (I^L,n+ RoiP,H)(t>R^rV 


( 66 ) 


Here, 




^ (sinh r) 


^—r^—ipr-\-inG 


^P±n 1 

V2 2 ’ 2 

for r —> oo 


ip — n ^ 1 

—-—;l + ip;-- 


sinh^ r 


(67) 


and 


pR,n 


= e“^(sinhr)-i+-^’('i-^,i 


ip — n . 1 


^—r ^ipr-\-in6 


for 


( 68 ) 


We see therefore that the normalisable eigenfunctions all vanish like (p (X e oc at large Meanwhile, 
eigenfunctions which are polynomial in T all are non-normalisable. Examples include 


ij=i 


cx (l + 2|vl/n, 


X (1 + 6|\E'P + el'I'l"^) 


(69) 


or 

oc (70) 

These ‘normal’ helds all correspond therefore to non-normalisable states, as can easily been checked since 

[ = ^. (71) 


To proceed, we now try to refine the connection between the monomer density and the non-compact 
degrees of freedom in the CFT. The operator coupled to AT —ATbn may be associated with a non-normalisable 
or with a normali sabl e state (for a discussion of the state-operator correspondence in CFT with non-compact 
targets, see, e.g., 61 ). In the first case, the conformal weight in the c = 2 -|- theory is negative, so the 


gap over the Ceff = z theory is 


C jU + 1) 1 2 

24 k-2 4(fc - 2) 24 ' 


(72) 


Such negative values of A would lead to exponents r < f, and there are no indications of such values in the 
numerical analysis. This means that the operator coupled to K — Kbn must be normalisable. While there is a 
continuum of choices (because of the continuous imaginary part of _)), it is reasonable, in the absence of other 
indications, to focus the bottoms of the continuum, corresponding to j = — |. Since monomer pieces involve 
moreover two links, the most conservative operators one can think of are associated to j = —1/2, n = 0 and 
j = —l/2,n = 2. From the exact form of the spectrum (beyond the mini superspace approximation (621), 


AlAl = + (n±fcw)^ 

^ k-2 Ak 


(73) 


we see in particular that operators with w ^ f) have too large a dimension to be relevant, and we do not 
consider them. 


22 














While the operator with j = —l/2,n = 0 has conformal weight zero with respect to the effective central 
charge, it is not necessarily trivial, as can be seen in the mini-superspace formulation, where it corresponds 
to a wave function F(l/2,1/2; 1, —|d>p). We propose therefore that the perturbation coupled to K — i^TeN 
can be represented by the action 

+ (K- if bn) [aF{l/2, 1/2; 1, -Id/p) + 3/2; 3, -Ivf-H] | (74) 

with some coefficients a,/?. For simplicity let us first look at the physics in the case where /3 = 0. 

/3 = 0 case. The hypergeometric function is monotonically decreasing from F = 1 for dt = 0 to F ~ 
when |'k| —>■ oo. If F < Fbn, we see that the action is minimised when F is at its maximum, so iTj = 0. This 
corresponds to a massive theory, with no critical degrees of freedom left. If F > Fbn, the action is minimised 
when F is as small as possible, so |4'| —>• oo. For F < Fbn, the classical action is kl « / d^x (F — Fbn)^ 
and for F > Fbn it is zero, so we expect a discontinuity of the first-order derivative, that is, a first-order 
phase transition. 

In the phase F > Fbn the fact that |d>| is infinite does not mean that the theory is trivial. In fact, this 
region is precisely the cylinder limit of the cigar, and one recovers the same free boson theory as the compact 
component of the critical point F = Fbn- It is reasonable to consider therefore the F > Fbn phase in 
our problem as a broken symmetry phase, and the compact boson describing the dense vertex model as a 
Goldstone boson. 

In fact, the phenomenology we propose is similar in spirit to the first-order phase transition occurring in 
systems described by a complex boson, with action 

A = J d^x + (F - FBN)|«'r(l - , (75) 

a model known to describe the superffuid transition (see, e.g., [62]). 

Generic case. More behaviours are possible when a and /3 are non-zero. Of course, the a term is always 
the most relevant and dominates whenever a ^ 0. When a = 0 exactly, on the other hand, the phe¬ 
nomenology is quite different from what we see in the F > Fbn perturbation. This is because the function 
a;^F(3/2,3/2; 3, —x^) has a maximum at a finite value of x. Irrespective of the sign of F — Fbn, the action 
is then minimised at this maximum, and for 20 = 0 or tt: both degrees of freedom then become massive. 
Since cos 20 is meanwhile relevant, it is very likely that this scenario persists in the presence of fluctuations, 
and does not lead to the observed first-order phase transition. 


We conclude from this discussion that the transition when F crosses Fbn in the untwisted model is 
driven by a field of dimension 0 with respect to the effective central charge Ceff = 2, and of dimension 
A = with respect to the central charge of the Black Hole theory. Unfortunately, the literature on 

perturbed field theories having non-normalisable ground states is lacunary, and we are not sure how these 
values of conformal weights are related with the scaling of physical quantities as F crosses Fbn- Certainly, 
both A = 0 and A = are very small in the domain of values of k reachable in regime III, leading to 

values of v very close to n = ^. The accuracy of our numerical calculations did not allow us to determine 
the possible difference ^ with reasonable accuracy. 

Turning now to the loop model, it is generally believed—and we have checked this numerically—that the 
ground-state energy of the theory (or the free energy per site in the 2D statistical point of view) is independent 
of the boundary conditions. This means that the scaling should be described by the operator with weight 


A = 


again, unless—as happens in the case of the six-vertex model and its loop formulation for 


4(fe-2) 

instance—the change of boundary conditions leads to the cancellation of some terms, and another, related 
but different, scaling dimension (see. 


e.g. 


63 


corresponding value 


for a discussion of this point). In our case however, the 

2(fc - 2) 


4fc-9 


(76) 
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seems to fit numerics reasonably well. Moreover, it is definitely the correct value for fc = 8,n = 0. In this 
case indeed, one can directly identify the operator coupled to K — with the two-leg operator, well-known 
in polymer theory, and whose dimension is A = ^ at this point. We conclude with the natural result that, 
to lowest order in the density, we have indeed that the monomer density is proportional to |5'p. 


6.5 The multicritical line 


The analysis of the perturbations from the Black Hole CFT at ©bn into the dense and massive phases can, at 
first sight, be extended to the perturbation along the multicritical line (the violet line in the phase diagram 
of figure]^. But while it is very natural to expect that this perturbation is relevant (and thus that there 
is a flow away from ©bn in all directions of the phase diagram in figure]^, the numerical evidence for this 
scenario is disappointing. 

Of course, investigating the critical behaviour along the multicritical line is a very complicated task. The 
precise location of the multicritical points away from ©bn and ©ds is unknown, and numerical estimations 
of critical exponents are plagued by the flow towards the dense or massive phases, as well as the proximity 
with the Ising or dilute critical surfaces. At this stage, we feel that an exploration of the phase diagram with 
conclusive numerical results is beyond the scope of this work. Here, we will content ourselves by investigating 
a much simpler question, namely what happens to the non-compact boson precisely at ©ds- This can be 
done by following continuously the levels (m,j) from ©bn- Moreover, we restrict for simplicity to the twisted 
case, where non-contractible loops have weight n. 

Before doing so, it is useful to recall a few results about the physics at the point ©ds- The latter belongs 
to a more general integrable line labeled ‘branch 0’ in 17 , and which was solved exactly by Bethe Ansatz 
Most of the results can actually be derived from the observation (reviewed above in section 
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that the configurations of the model with weights (n, ATdSjPdSi Hos) = |i0, 2) can be mapped onto the 

configurations of completely packed loops with weights n = n -|- 1, or equivalently, onto a Q = [n + 1)^ 
state Potts model. The latter model is critical for Q < 4, and well understood from a Coulomb gas 
construction 53 : the continuum limit is described by a compact boson, and parametrizing n = —2 cos ng 


the central charge reads 


c = 1-6-^ 
9 


(77) 


where eo is related to the weight Unc of non-contractible completely packed loops 
particular, for the ‘pure loop’ case where contractible and non-contractible loops 
eo = 1 — 5 so that 


by rinc = 2 cos ttbo . In 
have the same weight, 

(78) 


The point ©ds corresponds to n = 0, that is g = |, and so c = 0. Turning to the critical exponents, they 
read 


e(e-eo) gm'^ 

Xe,m- + 2 ’ 

(79) 

that is, at ©dsj 

e(3e — I) 

Xe,m — ^ + g - 

(80) 

We now come back to the levels (m,j). Contrarily to what was observed in the case of K perturbations, 
these appear to remain critical at ©ds- In particular, we have studied in detail the first few levels (0,j) in 
the sector of zero magnetisation. The following observations hold: 

• (0, 0) is still the ground state at ©ds (there is no level crossing); 

• (0,1) has at ©ds the conformal weight 

• (0, 2) has at ©ds the conformal weight A))'® = 2g. 
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From these observations it is quite natural to conjecture that the levels (0,j) become degenerate with the 
magnetic excitations of the compact boson at 0 dSj namely 


Kj = Xo,j 


(81) 


It is easy to imagine scenarios where the non-compact boson becomes massive, and the compact boson 
keeps its radius (as in the flow to the dense phase), or acquires a new value of the radius. It is also possible 
to imagine scenarios involving the emergence of Ising degrees of freedom, as observed on the blue surface 
in figure [2]— this would simply correspond, in the complex field picture, to the introduction of an additional 
gauge field [65| . But we cannot, for now, imagine a simple scenario connecting highly degenerate conformal 
weights at the 0ds point to the continuum of weights at the 0bn point. This will have to wait for further 
studies. Incidentally, this will require, in particular, furthering our understanding of the physics in the p = 0 
plane 
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7 Conclusion 


In this paper we have partly unraveled the mystery of the theta point in two dimensions. We have established 
that the BN point found in 17 is very peculiar from the point of view of critical systems: its correlation 


functions do not decay as discrete sums of power laws at large distance, but involve instead a continuum of 
critical exponents. This is expected to lead to extremely strong corrections to scaling, and explains why the 
results of numerical simulations have always been ‘off’ compared to the well-established theoretical values. 
We have also seen that the BN point is extremely unstable, and represents, in a certain sense, an infinite 
order of criticality. We note that these features are somewhat similar to the ones of the theta point for 


interacting trails. As discussed in 66 , this point is described also by a non-compact CFT, albeit of a rather 


trivial type (a Goldstone phase sigma model, free at low energy 67 ), with trivial critical exponents. 

One question that we have unfortunately not been able to fully answer is the role of the 0bn point in 
the phase diagram of polymers. We still believe that the point Oqs is the tricritical point for generic lattices 
and short distance interactions. It may be that it is also the generic tricritical point for the model we discuss 
in this paper, where loops osculate at the vertices (with weight r, see figure]^. This is what conventional 
thinking would suggest, since the critical theory at the 0bn point involves an infinity of relevant operators. 
Moreover, it is certainly difficult to imagine a flow from 0 ds— essentially described by a single boson—to 
0BN, which involves so many more degrees of freedom. On the other hand, we have not seen convincing 
numerical evidence for this flow. On top of this, it turns out that, for the model in figure the solvable 
point with universality class of 0ds enjoys an additional symmetry, which is broken as soon as ‘straight 
segments’ are allowed: the number of monomers on the even sub lattice minus the number of monomers on 
the odd sub lattice is indeed conserved by the interactions. This means, in fact, that the model enjoys an 
underlying SU{n) symmetry, larger than the 0(n) symmetry present in the rest of the phase diagram 36 


From this point of view, it is tempting to expect that there is an RG flow towards 0bn, which was suggested 
in the original paper 17 . But we have not seen clear numerical evidence for this either. It could be that 


there is a—yet unidentified—intermediate fixed point making these observations compatible. We hope to be 
able to report on this soon. 
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